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I. INTRODUCTION 

The fluctuation-dissipation (FD) theorem [lj-Q is a 
powerful tool in the study of equilibrium phenomena and 
has proved useful in a variety of circumstances. It re¬ 
lates two seemingly different phenomena in the system, 
namely, energy dissipation in the medium and the sta¬ 
tistical fluctuation of dynamical variables. The effect 
manifested itself in an unexpected manner in the under¬ 
standing of the Brownian motion of a particle as well 
as of the thermal noise in a conductor. In the case of 
the Brownian motion Einstein had shown that the coeffi¬ 
cient of diffusion was related to the coefficient of friction 
through a temperature dependent factor. Later Nyquist 
explained theoretically the experimental result of John¬ 
son that, in the absence of an applied current, the mean- 
square voltage in a conductor is related to the resistance 
of the conductor. Subsequently the FD theorem has been 
derived algebraically, starting from the first principles of 
equilibrium statistical mechanics (as well as in various 
other ways @^])- However, the symmetry principle or 
the conservation law, underlying such a powerful result, 
has not yet been fully clarified. Einstein had already 
noted qualitatively that one can understand the relation 
as follows. When a Brownian particle is subjected to a 
random force, the same force leads to two components - a 
statistical fluctuation and a drag - which must, therefore, 
be related. In this letter, we will show that the FD theo¬ 
rem can be understood as a consequence of the unitarity 
of the theory. Namely, we will derive the FD theorem 
starting from the unitarity of the S’-matrix. 

FD theorem holds in equilibrium, both classically as 
well as quantum mechanically. Here we will describe the 
quantum mechanical case setting h = 1 for simplicity. 
Let us recapitulate very briefly the algebraic derivation 
of the FD theorem within the context of a real scalar 
field theory. The derivation can be generalized to any 
other field theory in a straightforward manner. There 
are two essential components in the proof of the FD the¬ 
orem. The first (and the main) ingredient is that there is 
a relation between the correlated and the retarded prop¬ 
agators (Green’s functions) of the theory. This can be 
understood as follows. We note that, even though we can 
define two independent quadratic products from the ba¬ 
sic field operators, namely, (f(x°, x)^(0) and cj)(0)<j)fx o , x) 
(or, equivalently, the anti-commutator and the commu¬ 


tator), there exists only one independent thermal expec¬ 
tation value or ensemble average (this is also true at zero 
temperature for the vacuum expectation value). This is 
easily seen from the relation 

Tr (eP liH 4>(x°, x)0( 0)) = Tr {e~^ H (j){<S)(j){x 0 + iff, x)) . 

(1) 

This is known as the KMS condition [8j, [9j which follows 
from the cyclicity of the trace (as well as the identification 
that the Hamiltonian is the generator of time translation) 
and we have defined P = tpp- Therefore, if we define the 
correlated and retarded propagators as (G R (x) denotes 
the retarded Green’s function and stands for the 

ensemble average) 

C{x) = (f>(0)] + )p, iG R {x) = 9(x°)([<f>(x), 0(O)])/3, 

( 2 ) 

equation © leads, in momentum space, to the relation 
for the ensemble averages 

C[p) = -coth^ ImGjj(p). (3) 

We note here that the negative sign on the right hand side 
of © is a consequence of the field theoretic definition of 
G R (x ) in ©. The second element in the proof of the 
FD theorem is that, for weak external fields in a linear 
response theory, one can identify the response function 
x(p) with the retarded Green’s function G R {p) of the 
theory This leads to the result that the statistical 
fluctuations in a theory in equilibrium given by C(p) are 
related to the imaginary (dissipative) part of the response 
function through a temperature dependent factor. This 
relation is known as the FD theorem. 

We note for future use that the time ordered product 
of two field operators which defines the Feynman propa¬ 
gator of the theory can be written as 

T(<j>(x)<j>( 0)) = 6(x°)cp(x)cp( 0) + 0{-x°)(j>(Q)cj){x) 

(4) 

which, upon taking the ensemble average, leads to 

iG f (x) = C(x) + l - ( G r {x ) + G r (-x)) . (5) 
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This ties in with the Einstein observation within the field 
theory context, namely, C{x) and iG R (x) correspond to 
two parts of the same Feynman propagator and the FD 
theorem (J31) only gives a relation between these two com¬ 
ponents. However, it does not explain why this special 
relation should actually hold. We also note from f4]) that 
since the anti-commutator is Hermitian while the com¬ 
mutator is anti-Hermitian, it follows that 


and the two sets of matrices in m and m are related 
through a 2 x 2 unitary matrix (see, for example, eqs. 
(2.42)-(2.43) in 13). At finite temperature, the results 
of chapter 5 in JT^l show that a cutting description holds, 
say for example, as a matrix for the self-energy, namely, 
element by element so that we can write 

2ImS ab (p) = U^ b (p) + Yl* b (p), a,b = ±, (11) 


Im G f = —C. (6) 

This relation holds in coordinate as well as in momentum 
spaces since both ImGf(a;) and C(x ) are real and are 
even functions of x. As a result, we can also write the 
FD theorem ([3]) as 

Im Gf(p) = -C(p ) = coth lmG R (p). (7) 

We will use this form of the FD theorem and show that 
the origin of this powerful relation lies in the unitarity of 
the theory. 


II. UNITARITY AND THE FD THEOREM 


In a quantum field theory, conservation of probability 
is encoded in the unitarity of the 5'-matrix, 


S'S = SS^ = 1 , ( 8 ) 


which implies that the imaginary part of any amplitude 
in the theory can be given a cutting description. At zero 
temperature, such a cutting description, even though not 
required by unitarity, holds graph by graph and is known 
as the Cutkosky rule |10|. At finite temperature, it has 
also been shown 0 [l2| to all orders within the closed 
time path formalism, that such a description indeed holds 
for the imaginary part of any amplitude as a whole (and 
not graph by graph). 

The closed time path formalism 00 , which we will 
use, is ideal for the study of nonequilibrium phenomena. 
It involves doubling the field degrees of freedom (thereby 
doubling the interaction vertices as well). (We follow the 
notations and conventions of fl2j | and make use of many 
results derived in chapter 5 there in order not to dupli¬ 
cate technicalities.) As a result, the Feynman propagator 
(Green’s function) and the self-energy become 2x2 matri¬ 
ces labelled as G a b, F 0 &, a, b = ± (± denote the doubled 
thermal degrees of freedom) of the forms (see chapter 2 
in 0) 



We note here that the Feynman Green’s function Gf 
which we have already defined in ([5]) coincides with G++. 
One can also define 2x2 matrices incorporating physical 
propagators (Green’s functions) and self-energies of the 
forms 


G = 


( 0 
\G r 





( 10 ) 


where n ^ b (p) and n ^ b (p) represent the cut diagrams for 
the self-energy with the cut towards the left and right 
respectively as shown in Fig. |T] The retarded self energy 
is given as the sum 

Z R (p) = S++(p) + E+-(p). (12) 

So, for example, for the imaginary parts of the Feynman 
and the retarded self-energies we can write 

2ImS ++ (p) = n^ + (p) + n£ + (p), 

2imE fl (p) = n^ + (p) + n^_(p) + n^ + (p) + n£_(p) 
= n£ + (p) + n£_(p), (13) 

where we have used the result from 0,0| that the right 
handed cut graphs cancel in the retarded self-energy. To 



FIG. 1: Left and right cut self-energy diagrams. 

proceed further with the calculation of the imaginary 
parts of the self-energies in (fl3l) . we collect various prop¬ 
erties satisfied by the cut diagrams n % b (p) and n^ b (p). 
First, we note from Fig. 1 that 

U ab(~p) = Hf a (p), (14) 

which, in fact, holds graph by graph, but (fT4l) suffices for 
our purpose. It is also known that the sum of all the cut 
diagrams are real, namely, 

(h^(p))* = n L ab ( P ), (n« (p))* = n* (p). (15) 

Furthermore, there is an underlying symmetry in individ¬ 
ual graphs for self-energy at finite temperature, namely, 
under (a, b) (—a,—b),p —>• — p and complex conjuga¬ 
tion, any self-energy graph remains invariant (there is a 
corresponding underlying symmetry for any graph, not 
necessarily the self-energy, in a thermal field theory). 
Together with m , this leads to the result n ^ b {p) = 

( n -6,- a (p)) . Equation (flbl) then leads to 
u ab(p) =n\- a (p)i a,b = ±. 


(16) 
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We also use a crucial theorem (proved in chapter 5 of 
Pj ) that the sum over the thermal indices on the cut 
side of the self-energy diagrams vanishes 

E n a6W=0 = E n a6(p), (17) 

a=± 6=± 

which leads to 

n* + (p) = -n?_(p) = -n£_(p), (is) 

where we have used m in the last step. Using the rela¬ 
tion GBD , we can now write (113l) in the simpler form 

2 imS ++ (p) = n l + ( P ) + n ++{p), 

2ImE fl (p) = n L ++ {p) - n« + (p). (19) 

To relate ImE_|_ + (p) and ImEij(p), we need one final 
ingredient. Let us consider the closed loop graph shown 
in Fig. [2]with a II ^b(p) insertion (a, b fixed) 


V 



FIG. 2: Closed loop with a Yl^ b (p) insertion. 


lab = J ^n«(p)G a ,_ a (-p), ( 20 ) 

where we have used the property of thermal graphs (see 
eq. (5.41) in pj l that when one end of a propagator is 
in the cut side of a graph, it is completely determined by 
the thermal index of the uncut vertex. Furthermore, if 
we let p —> — p in the integral and use m, this leads to 


which, in particular, implies that 

n?+(p) = e-^n L ++ {p). (24) 

This equation ensures that, in the zero temperature limit, 
the self-energy diagram with a right cut vanishes when 
Po > 0, as is known from the study of cutting rules at 
zero temperature. Substituting equation (1241) into m 
we obtain a direct relation between the imaginary parts 
of the Feynman self-energy and the retarded one, 

ImE ++ (p) = coth^ ImSjj^). (25) 

This relation follows from the unitarity of the theory and 
has the same form as 0, but holds for the self-energies. 

To make connection with the FD theorem, let us recall 
that the Green’s function and the self-energy in 0 are 
related as 


G- 1 (p)={g^(p)) '-ECp), (26) 


where G'°> denotes the free (tree level) Green’s function 
of the theory. The 2x2 matrix propagator in © can 
be simply inverted (in the space of thermal indices) and, 
with some analysis, (EU1) leads to the relation 


ImE ++ (p) = - 


det G(p) 


Im G ++ (p) 


+ 


1 


det G(°)(p) 


Im g£]_(p). 


(27) 


In a completely parallel manner, one can determine from 
m the relation 


ImE R (p) = - 


det G{p) 


Im Gr{p)+ 


det G(°)(p) 


ImG^V). 

(28) 


From the fact that G(p) and G(p) (as well as G^(p) 
and G(°) (p)) are related by a unitary matrix (so that 
their determinants are the same) and using (l25l) together 
with relation (2.51) in P] we obtain the FD theorem 0, 
namely, 


lab = J d 4 pn£ a (p)G a - a {p) 

= j d 4 pIlt(p)e- afSp °G a> - a (-p). (21) 

Here we have used the fact that when a thermal prop¬ 
agator has two opposite thermal indices, it satisfies (see 
eq. (5.43) in 0) 

G a - a {p) = e~ a ^ Po G a - a {—p), o = ±, (22) 

which follows from the KMS condition. Comparing EU1) 
and CD we obtain the relation 

nf b (p)=e- Q ^n^ a (p), a,b = ±, (23) 


Im G F (p) = lmG ++ {p) = coth ^- Im G R (p). (29) 

This gives a derivation of the FD theorem starting from 
the unitarity of the theory. 

III. CONCLUSION 

Our derivation shows that the physical content of the 
fluctuation-dissipation theorem may be understood as 
arising from the unitarity of the 5'-matrix. As we have 
noted earlier (see 0 as well as the discussion following 
0), this result expresses a general relation between the 
fluctuating properties of a system in thermal equilibrium 
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and the response of the system to a weak external per¬ 
turbation. Finally, we would like to comment on the clas¬ 
sical limit of our result. In our entire analysis, we have 
set ft = 1. Restoring the factors of ft, the temperature 
dependent factor in (l29l) takes the form 

coth ^ ^ ft coth . (30) 

2 2 

In the classical limit, ft —>• 0, this factor reduces to ■ 3 =- 
. . ... ftpo 

which leads to the classical fluctuation-dissipation rela¬ 


tion in © or HMD. 
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